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Abstract. This paper studies theoretically the ultraslow light phenomenon in 
BoseEinstein condensates of alkali-metal atoms. The description is based on the linear 
approach that is developed in the framework of the Green function formalism. It is 
pointed out that the group velocity of light pulses that are tuned up close to resonant 
lines of the alkali atoms' spectrum may strongly depend on the intensity of the external 
static magnetic field. The possibility of filtering optical pulses using the ultraslow light 
phenomenon in Bose condensates is discussed. 
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1. Introduction 

Bose-Einstein condensate (BEC) is one of the most impressive examples when the 
matter demonstrates its quantum nature on the macroscopic level. Now this system is 
interesting also due to the possibility to observe the electromagnetic pulses propagating 
with extremely slow group velocities in it [1] . Evidently, the choice of alkali- metal vapors 
corresponds to the fact that a behavior of a many-particle system consisting of alkali 
metals is more predictable and convenient from the standpoint of quantum mechanics 
than a behavior of the system consisting of atoms with more complicated structure. The 
main reason for that is a hydrogenlike structure of alkali metals. Therefore, the energy 
spectrum of these atoms and interaction peculiarities with an external electromagnetic 
field are well-studied to the present moment. 

For the many-particle system consisting of alkali atoms one can construct rather 
exact microscopic description, which includes the effects of collective interaction of 
these systems with electromagnetic waves. This description is based on an approximate 
formulation of the second quantization method [2] in the case of a presence of bound 
states of particles. In the framework of this formulation one may consider alkali metals 
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as bound states of two "elementary" particles: the valence electron and the atomic core. 
This method allows not to consider an atom as a compound particle, but describe it as 
an elementary object with preserving the information relating to its internal degrees of 
freedom (wave functions, energy spectrum, etc). 

On the basis of the developed microscopic approach in Refs. j3j HJ E] it was studied 
the response of a system to an external electromagnetic perturbation in the framework 
of the Green-function formalism. In terms of the Green functions the expressions for 
the permittivity and magnetic permeability of an ideal gas of hydrogenlike atoms with 
a Bose condensate. It was found that the permittivity demonstrates the resonance 
behavior at frequencies close to the allowed transitions between certain quantum states 
of an atom. Note that this behavior of the permittivity is well-known from the classical 
approaches (6j. But, the proposed theory allows to estimate the values of all microscopic 
quantities that appear in the expressions characterizing the response of a gas with a 
BEG 

It it easy to see that the energy spacing between different states of an atom can 
strongly depend on the intensity of an external magnetic field. Aside the cases, in 
which this field can be specially turned on, it can appear due to the usage of magneto- 
optical traps and bias fields in BEC-related experiments. The mentioned dependence 
of the energy spectrum on the magnetic field corresponds to the Zeeman splitting (or 
the Paschen-Back splitting) of the ground and excited states of an atom. In turn, the 
resonance behavior of the refractive index is one of the necessary conditions for the 
pulses slowing in a Bose condensate jl]. Thus, one can expect that by the change of 
the intensity of a bias field it is possible to control the propagation velocity and to filter 
light pulses that propagate through a BEC cloud. 

2. Linear response of a gas to an external probe radiation 

Firstly, it should be noted that in the present paper the description is based on the weak- 
field approximation. In other words, it is assumed that an external electromagnetic field 
(laser pulse) is probe and practically do not influence on the occupation of the states in 
the system. This limitation is necessary for the usage of the results of the linear-response 
theory, which is based on the Green-functions formalism [7]. The mentioned fact 
excludes the possibility to describe the influence of an additional coupling laser (quantum 
interference effects) on the system under consideration. Note that the coupling field 
was used in the experiment pLJ to provide an electromagnetically-induced transparency 
(EIT) at the resonance regions. But, as it is shown in Ref. jl], the usage of quantum 
interference effects is not a necessary condition for the pulses slowing. Moreover, the 
coupling field may influence on the energy spectrum of the atoms. Thus, it inevitably 
produce additional difficulties in the construction of the description. 

Usually, for the achievement of the BEC regime in the systems consisting of alkali- 
metal atoms, the dilute vapors with the density of particles v < 10 14 cm -3 are used. 
These vapors with a good precision can be considered as ideal gases [H]. Hence, to 
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describe microscopically the propagation of electromagnetic pulses in this medium, it 
is convenient to use the results [3]. There it was shown that the expression for the 
permittivity of the gas in the BEC state can be written as follows: 



a 



e-^k, u) « 1 + ^ V ^ fl 7 JWWI , (1) 
y ' ^ u + Ae a0 + i laP v ; 

where w, k are the frequency and wave vector of the external field (probe laser), v a is 
the density of condensed atoms in the a quantum state, Ae a p is the energy interval 
between the a and /3 states, and 7 a/ g is the linewidth related to the probability of a 
spontaneous transition between these states. The quantity er a/ g(k) is the matrix element 
of the charge density, which is defined by the relation 

/r / m \ / m \ ' 

rf y<^(y)<My) [ ex P y— k y) - ex P ( _i ~ ky JJ > ( 2 ) 

where ip a {y) is the atomic wave function in the a state, e is the electron charge 
absolute value, m p and m e are the masses of the atomic core and electron, respectively 
(m = irip + m e ). In the case of the allowed dipole transitions, in the linear order over 
(ky) <C 1, one gets 

W(k) « ikd a p, d a(3 = ej ydy<p* a (y)<pp(y), (3) 

where d a p denotes the atomic dipole moment, which is related to the transition a <H- 0. 

Considering that the density of the excited atoms is negligibly small in comparison 
with the density of atoms in the ground state (the case of a low intensity of the probe 
laser), one can neglect of the second term in Eq. ([1]). Therefore, one gets the expression 

e -1 (k, w; B) ps 1 + ^ a Q/3 [u + Ae a p(B) + i^]' 1 . (4) 

Here the set of quantum numbers of ground and excited state are denoted by the a and 
13 indexes, respectively. As in this paper the limit of zero temperatures is considered, 
the quantity 7 a/ 3 corresponds to the probability of a spontaneous transition between the 
/3 and a states. Note that by the term "ground state of an atom" relates to all sublevels 
of the ground state (for alkali metals these levels have zero orbital number), which are 
splitted by the hyperfine and Zeeman interaction. According to Eqs. (0Q) and ([3]), the 
quantity in Eq. (Tjfl) 

a aj3 = iiru a d 2 a p/3 (5) 

defines the dependence of dispersion characteristics on the atomic polarization and 
the density of atoms in the BEC state. Note also that the dependence of dispersion 
characteristics on the intensity of a static magnetic field B (see Eq. (jJJ) is mainly defined 
by the energy interval Ae a p{B) between corresponding ground and excited states of an 
atom. Thus it is necessary to study dependencies of the resonant transitions on the field 
intensity in detail. 
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3. Hyperfine structure of the sodium D 2 line in the external magnetic field 

Now it is important to obtain the energy spectrum of alkali-metal atoms in the presence 
of a magnetic field with account of the hyperfine interaction. To this end, one should 
consider the following Hamiltonian (see also Ref. [9]): 

V = A hfs I-3 + B hf8 Q + ^ B {gilz + gjJ z )B, (6) 

where A^fs and -B/i/s are the hyperfine constants related to the magnetic dipole 
and electric quadrupole interaction of the nucleus and valence electron, respectively 
(Bhfs = for the states J = 1/2). The quantities I and J are the operators of the 
nucleus and electron total angular momentum, respectively, I z and J z are the operators 
of the projections of these operators on the magnetic field B direction, \lb is the Bohr 
magneton (in the units h = 1), and gj and gj are corresponding Lande factors. The 
operator Q can be expressed in terms of the operators I and J as follows: 

a = 3(i-j) 2 + fi-j-/(/+i)j(j + i) 

W 2/(2/ - 1)2 J(2J - 1) ' 1 ' 

Note that for the splitted levels of the ground state the spectrum can be found 
analytically. At the same time, for the excited levels (triple and more degenerated 
states in the total momentum projection mp, F = I + J) it is rather difficult to find 
corresponding analytical expressions. But in this case one can find it numerically. One 
may verify that the limiting cases of the numerical results lie in a good agreement with 
the analytical formulas. 

By the example of a sodium atom, it is possible to study the dependence of D 2 - 
line components on the magnetic field in detail. This line corresponds to the transition 
between the 3<Si/2-state sub levels and 3P3/2-state sublevels. It is not difficult to obtain 
the Zeeman splitting of these levels with account of the hyperfine interaction (see Eqs. ([6]) 
and (J7|)) as it is done in Ref. [9]. As one can see, the resonant D 2 line splits to a large 
number of components. Each component can be found from the difference between 
corresponding excited and ground state with the use of the selection rule for the quantum 
numbers F and mp (AF = 0, ±1; Amp = 0, ±1). 

Usually, depending on Amp value, one can say about linearly or it -polarized 
transitions (Amp = 0), and circularly or a ± -polarized transitions, Amp — ±1 [6]. 
Note that in the experiments with ultraslow light the laser with the certain polarization 
is used. Thus, one needs to consider the transitions between the states of an atom 
corresponding to the certain polarization. Hence, in the next calculations the response of 
the system to the linearly-polarized pulse is considered. For this radiation the resonance 
frequencies correspond to the transitions AF = 0, ±1; Amp = 0). The magnetic-field 
dependence for these components can be found from the dependencies for the Zeeman 
splitting of the ground and excited states (see, e.g., Ref. [5]). Accounting the mentioned 
selection rules, one can find the dependencies for the "upper" ^-transitions (transitions 
from the ground state with F = 2) that are shown in Fig. [TJ 
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Figure 1. Components of the sodium Di line in the external static magnetic field 
(F = 2; AF = 0,±1; Amp = 0). The number in brackets corresponds to the 
momentum projection mj? of the quantum state. 



4. Dependence of the group velocity on the magnetic field intensity 

Now one can conclude that at certain intensities of magnetic fields at the same frequency 
two different transitions are possible (regions of the lines intersection, see Fig. [p. These 
regions may be very interesting from the standpoint of the experiments related to 
the ultraslow-light phenomenon. The reason is that the value of the group velocity 
depends on the dispersion characteristics of the medium. In particular, it depends 
on the derivative of the refractive index by the frequency. Note that the smaller 
the difference between two resonant frequencies, the narrower the spacing between 
corresponding peaks of the refractive index, and the steepness of the slope is larger. 
But this characteristic does not mean that in the region of intersection one should 
expect the maximal reduction of the group velocity. In this case, the natural linewidth 
of the excited state component (sublevel) plays an essential role. These quantities for 
the excited sublevels can be obtained from the relative hyperfme transition strength 
factors Sff' [9], Iff 1 = SW'F e /2, where for sodium atom T e = 9.795 MHz. For 
example, if the atoms are prepared in the ground states with F = 2, one must take 
corresponding strength factors for the "upper" transitions (see Fig. [T]), S 2 i = 1/20, 
S 22 = 1/4, and S 23 = 7/10. 

Next, analyzing dependencies that are shown in Fig.[U one can choose the frequency 
of the linear-polarized light close to the lowest intersection of the eigenfrequencies 
u)\ (dotted line in the lower part denoted by (0)) and u 2 (dotted line in the lower 
part denoted by (-2)). In particular, it is convenient to take co^ = coo + A, where 
Uq is D 2 -line resonance frequency without account of the hyperfme interaction (oj = 
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508.848716 THz), and A is the chosen detuning (A = —45.0 MHz). It is also necessary 
to consider that two corresponding ground-state sublevels are occupied (quantum states 
|1) = \F = 2,mp = 0) and |2) = \F = 2,mp = —2)). In this case one can simplify 
calculations by accounting only two resonant summands in (T4]). Really, as one can 
conclude from Fig. (TJ in this case the major role play the transitions corresponding to 
the eigenfrequencies Ui and uji- Therefore, one gets rather simple expression for the 
permittivity, 

2 

e-\k,uj;B) ^l + ^ aj [uj -uo J {B)+i lj }- 1 . (8) 
i=i 

where the quantities jj = SFF'F e /2, aj = 4iih'jd FF ,/3 (see Eq. ((5])), d 2 FF , ~ 
Sj?i7v(3.52er ) 2 [9], r is the Bohr radius, and e is the elementary charge. The summand 
j — 1 in this case corresponds to the transition quantum numbers F — 2, F' — 1 and 
the summand j = 2 corresponds to F — 2, F' — 2. 

By introducing the real and imaginary parts of the permittivity (e' and e", 
respectively) it is easy to find the relations defining the refractive index and damping 
factor quantities: e' = n' 2 — n" 2 and e" = 2n'n" . In the transmission region the refractive 
index n' defines the group velocity of the signal, and the damping factor n" defines the 
intensity of the transmitted light, 

v g (u, B)=c {n\u, B) + w[dri(w, B)/du}}~ 1 , (9) 

B) = I (u) exp [-n"(w, B)kL), (10) 

where Io{&) is the spectral density of the initial pulse, L is the characteristic linear size 
of a condensed cloud. 

Next, by setting u = u^ e \ v\ = V2 = 2 x 10 10 cm -3 , and L = 0.004 cm, it is easy 
to calculate the quantities ([9]) and ( fTDl for the different values of the intensity of a bias 
field. By the use of these results one can get the dependencies that are shown in Fig. |2j 
As it is easy to see, these dependencies have a rather nontrivial form. One can conclude 
that in some regions the velocity can be reduced more than in ten times by the change 
of the field intensity only in several gauss. Moreover, in some regions (\u^ — Ui\ < ji) 
the group velocity has negative values. But, as it is shown in the lower part of Fig. |2l 
in these regions the signal is mostly absorbed. The absorption rate can be made much 
lower by the use of a gas with the lower density or with smaller characteristic size L. 
In this case one should observe the light pulses propagating in a BEC with a negative 
group velocity. Note also that the propagation of the light pulses with a negative group 
velocity for similar systems is well-studied both theoretically and experimentally (see 
Ref. [10] for details). 

5. Possibility to filter the signals by the use of the ultraslow-light 
phenomenon in a BEC 

To show a possibility to filter optical pulses in a BEC in the presence of the external 
static magnetic field it is convenient to choose some model. In particular, when the 
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Figure 2. Dependencies of the group velocity and relative intensity on the external 
static magnetic field for the resonant light pulses with the constant detuning A = 
-45.0 MHz. 

spectral density of the initial pulse has a normal (Gaussian) distribution, it can be 
written as follows: 

J (w) = -A= exp [-(u - u ) 2 /2a 2 } . (11) 
a\/ZTT 

To get appreciable results, it is convenient to choose u = 508.848052 THz, a = 50 MHz 
(the signal can effectively interact with several hyperfme components of the sodium D 2 
line, see Fig. [1]), and the linear polarization of the signal. 

Next, it is necessary to choose the gas properties. For example, one can take the 
gas of sodium atoms with two occupied ground states: |1) — \F — 2,mf = 1) and 
|2) = \F = 2,mp = 2). The resonant transitions dependencies on the magnetic field for 
these states are marked by solid bold lines in Fig. [TJ Next, taking for the convenience 
Ui — v-i — 10 12 cm~ 3 , B = 25 G, and L = 40 /zm, in accordance with Eqs. (jSJ) and ( fTUl) . 
one can get the dependencies that are shown in Fig. [31 

Therefore, from Fig. [3] one can see that different parts of the pulse have different 
time delays resulting from the fact that the group velocity of these parts in gas a with 
condensates has different values due to the ultraslow-light phenomenon. It should be 
noted that that at certain frequencies the signal is mostly absorbed by medium. But 
this is not enough (see dashed line in the lower part of Fig. [3]) to reach good parameters 
for filtering optical pulses. One needs also to use additional techniques for a registration 
of the signal only in the certain time interval. As it easy to see from Fig. [3] (bold solid 
line in the lower part corresponds to the time interval from 12 to 14 ns), by the use of 
the mentioned method one can get filtered pulses with relatively good parameters. 

In conclusion, it is important to note that the cases, which is studied in this paper, is 
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Figure 3. Frequency dependencies of the time delay and relative intensity of the signal 
that propagates through a condensed gas with the occupied states \F = 2 : mp = 1) 
and \F = 2,m F = 2). 



demonstrative in some sense. In the framework of the proposed approach one can obtain 
the results for another polarizations, occupied states, another frequency detunings, and 
also another types of alkali atoms. Evidently, this fact expands the possibility for the 
experimental observation of the effects that are studied in this work. 
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